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Abstract 

Wc analyze a simple asset transfer model in which the transfer amount is a fixed fraction / of the giver's 
wealth. The model is analyzed in a new way by Laplace transforming the master equation, solving it 
analytically and numerically for the steady-state distribution, and exploring the solutions for various values 
of / G (0, 1). The Laplace transform analysis is superior to agent-based simulations as it does not depend 
on the number of agents, enabling us to study entropy and inequality in regimes that are costly to address 
with simulations. We demonstrate that Boltzmann entropy is not a suitable (e.g. non-monotonic) measure 
of disorder in a multiplicative asset transfer system and suggest an asymmetric stochastic process that is 
equivalent to the asset transfer model. 
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1. Introduction 

A vibrant research theme in econophysics is the analysis of asset exchange models. In these models, a 
large number of agents iteratively exchange assets, typically representing monetary amounts. In the simplest 
model that has been considered, the transfer amount is constant and independent of the agent's wealth, 
producing an exponential wealth distribution in the steady state (see [l[ for a review). More complicated 
fractional exchange models have also been considered by several authors [D, Q , in which the size of each 
transfer is a linear function of the wealths of the agents involved in the exchange. 

It has been found both analytically and numerically that the steady-state wealth probability distribution 
function pJw) in fractional exchange models depends strongly on the parameters that characterize the 
exchange [J]. Certain parameter values or exchange rules yield a strongly peaked distribution with an 
exponential tail, while other values yield a broad distribution with Pareto-likc qualities. The dichotomy is 
exemplified by two simple models. If the transfer amount is a fixed fraction / of the giver's wealth (the 
giver is the agent who surrenders the asset in the transfer), then the resulting steady-state distribution is 
strongly peaked and decays exponentially in the tail. If, on the other hand, the transfer amount is a fixed 
fraction / of the poorer agent's wealth, then one finds a broad steady-state distribution, which can be fitted 
well by a power law with exponent —1 across a broad interval of wealths. In this paper, we refer to these 
two models as the giver scheme and the poorer scheme respectively!]] 

In some of the asset exchange models considered in Q and several other studies Q, Q , the fractional 
exchange amount is a random linear combination of the wealths of the participating agents. The controlling 
parameter is the saving propensity, A, which determines the fraction of the agents' wealths that they do 
not offer to exchange. Comparing the output of simulations for the giver scheme and the exchange schemes 
based on the saving propensity, one observes that the schemes are closely related, with / ~ corresponding 
to A w 1. If the saving propensity is the same for all the agents, then the resulting steady-state distributions 
are similar to those obtained for the giver scheme. On the other hand, if the saving propensity is uniformly 
distributed, the steady-state distribution is a power law, p s (w) oc w~ 2 . In a recent study based on numerical 
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simulations [7[ , it was found that a combination of the poorer and giver schemes in one simulation results in 
a power-law wealth distribution whose exponent depends on the relative contributions of the two schemes. 
The more agents follow the giver scheme, the greater the exponent. 

Asset exchange models can be treated analytically via a kinetic or master equation, which tracks the 
rate of change of the number of agents at any given wealth. In particular, the master equation for the 
giver scheme has been derived in [8|. These authors found an expression for the second moment in the 
steady state, which agrees with the expression found in for a similar model by assuming the gamma 
distribution of wealth. The standard deviation converges to its steady-state value exponentially on a time 
scale ~ [f(l — /)] • The authors also found the asymptotic behaviour of the wealth distribution at small 
values of wealth. The master equation for the poorer scheme was derived recently [loj ]. but its solutions have 
not yet been studied. In [f|, the authors derived the kinetic equation for the case of uniformly distributed 
saving propensity. They demonstrated that the solution follows a power law with the same exponent as in 
the simulations. The kinetic equation approach was also used in to analyze self-similar solutions of a 
non-conservative asset exchange system. The author found a closed-form solution in the limit of continuous 
trading by means of the Laplace transform and observed that the distribution exhibits power-law behaviour 
at large wealths. 

The dependence of the relaxation time on the exchange parameters has been investigated numerically 
in [l2| for the models considered in 0, [fjj. The relaxation time-scale was found numerically to scale as 
~ (1 — A) -1 . This is consistent with the values found analytically in the giver scheme for the standard 
deviation. The authors also considered how the relaxation time depends on the number of agents but failed 
to find any significant trend. 

Recently, much effort has been directed profitably at developing more sophisticated and realistic multi- 
agent models to be analyzed by means of numerical simulations. In the present paper we take the opposite 
tack and return instead to one of the simplest multiplicative models, the giver scheme. We show that its 
master equation can be solved efficiently by a Laplace transform technique. Armed with this new tool backed 
by multi-agent simulations, we identify the following new properties of the system. (1) We get precise values 
of various quantities such as the steady-state entropy as a function of the model parameter /, independently 
of the number of agents. (2) We explore the thinly studied regime 1/2 < / < 1 and identify its unusual 
properties, e.g. oscillations in p s (w). (3) Using multi-agent simulations, we investigate how the Boltzmann 
entropy evolves with time as the system approaches equilibrium and argue that the Boltzmann entropy is 
not a suitable entropy for the giver scheme, even though the system is closed and conservative. (4) We 
propose a simple asymmetric stochastic process that is equivalent to the giver scheme. (5) We investigate 
how the degree of inequality, characterized by the Gini coefficient, depends on /. (6) Finally, we apply 
phase-space techniques from statistical mechanics to the giver scheme in order to illuminate the difficulties 
and opportunities that this asset transfer model presents. 



2. Giver scheme 

We consider a simple asset transfer model, in which the transfer amount is equal to a fixed fraction of 
the giver's wealth^ If w g is the giver's wealth and w r is the receiver's wealth prior to the transfer, then 
their wealths after the transfer are given by w g — Aw and w r + Aw respectively, with Aw = fw g and 
/ G (0, 1). The model comprises a large number of agents, who are assigned wealths initially according to 
some distribution. The transfers are assumed to take place over a fixed time interval At. At each discrete 
time ti, the agents are divided randomly into pairs and the transfer formula is applied to each pair. The 
transfers are complete by the time t,+i = U + At and the process repeats at the time U+i. The probability 
of drawing any pair is the same. In each pair, the giver is assigned randomly regardless of the wealths of 
the agents. 

The master equation for this system was derived in Q and is given by 

dp(w,t) , . 1 ( w \ 1 f w , . fw-w 1 \ . . , 

- -p(w,t) + — -pi- -t)+— dw' p[——,t)p(w',t). (1) 



dt ^ ' 2(1- fy vi-/' ) If h \ f 



2 The giver is also called the payer or the loser in the literature. 
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It is easy to verify that the mean of the distribution, [i\ — dw wp(w, t), docs not depend on time. Upon 
integrating by parts, one arrives at the evolution equation 



dt 



= -/(I - j>2 + M (2) 



for the second moment, H2(t) = / °° dw w 2 p(w,t), first reported in Q. This equation can be solved for the 
variance 

a\t) = M2 (t) - nj = ( M2 (0) - ^L) e-'t 1 "')* + (3) 

In the steady state, one has a s = a(t —> oo) = — Z)] 1 / 2 . For simplicity, we assume henceforth that 

the mean of p(w, t) equals unity|f| In the following sections, we are mostly concerned with the steady-state 
distribution p s {w) = p(w, t — > oo). 



3. Laplace transform of the master equation 

The Laplace transform of the master equation in the steady state is given by 

g{z) = \g{z-fz)+ l -g{z)g{fz), (4) 

with g(z) = dw e~ zw p s (w). Note that the functional equation (0]) applies to any integral transform 
whose kernel depends only on the product of the arguments of the function and its transform. For / = 1/2, 
the functional equation has a closed form solution 

where C is a complex-valued constant. Using the definition of the transform, we have <?(0) = 1 from the 
normalization of p s (w) and <7'(0) = — 1 from the assumption that p s (w) has unit mean, which gives (7 = 1. 
Applying the inverse Laplace transform to this solution gives the exponential distribution, p s (w) — e~ w , 
which was obtained in [8| by substituting simple "test" functions into the master equation. No closed-form 
solutions have been found for other values of / E (0, 1). 

The Taylor expansion of g(z) at z = can be derived by substituting the expansion in Q and using 
g(0) = 1 and <?'(0) = — 1. For the first four terms of the expansion, this procedure gives 

In general, for g(z) — X^^Lo a n(—z) n /n\, clq = 1, and a% = 1, we obtain 

a n = y( n ) fk<lka r k ; fom>l. (7) 

few 1 - /" - (i - /)" 

Since g{—z) is the moment-generating function for the distribution p s (w), the n-th moment of the distribu- 
tion /j, n equal a n , i.e. all moments of the steady-state wealth distribution can be computed for any / using 
the recursive formula ([7]). 

Using the Taylor expansion, one has a n — > 1 as / — > and hence g(z) —> e~ z . Note that the functional 
equation Q becomes an identity for / = and g(0) = 1. Taking the inverse Laplace transform of g(z) = e~ z , 
formally one gets p s (w) = 5(w — l). However, this wealth distribution is never reached because the relaxation 
time scale t r = [f(l — J)] -1 determined from tends to infinity in this limit. Indeed, p s (w) is equal to 



If the mean m 1, one can consider the function q(x) = fj,ip(fj,ix), which is normalized and has unit mean. 
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(a) (b) 

Figure 1: The Laplace transform g(z) for / = 0.1 obtained by solving ||4} iteratively. |(a)| Rc[g(z)] (top left panel), Im[g(2:)] 
(bottom left panel), and |g(z)| (right panel) versus r along the real axis (dashed curve), the imaginary axis (solid curve), and 
the line inclined at 9 = 60° to the real axis (dotted curve). The variables r and 9 are defined by z = re lB . |(b)| A view of the 
real (top) and imaginary (bottom) parts of g(z) (values above 1 and below —1 have been cut off). 



the initial distribution if / = 0. On the other hand, the Taylor expansion does not have a limit as / — > 1. 
The functional equation Q has the solution g(z) = 1 when / = 1, but it does not satisfy the condition 
g'(0) = — 1. It appears that p s {w) does not have a proper limit as / — > 1. Note also that the relaxation time 
tends to infinity as / approaches unity as well. 

The asymptotic behaviour of g(z) at infinity can also be deduced readily from the functional equation. 
In the directions in the complex plane for which one has g(z) — t as \z\ — > oo, the equation 



j(z - fz) = 2g{z) 



(8) 



must be approximately true for large enough \z\. Assuming a power-law shape \g(z)\ oc \z\ a as \z\ — > oo, 
equation © gives 

<*= -, : x (9) 



log 2 (l-/)' 

By Watson's lemma [l3j . this is consistent with the asymptotic behaviour p(w) oc w' 
found in by the method of dominant balance. 

The functional equation Q can be solved iteratively when it is cast in the form 



9i(z ~ fz) 



as w — > that was 



(10) 



where g%(z) is the i-th iteration. Experimentation shows that the choice go(z) = 1/(1 + z) works well for all 
/. A detailed description of the computational procedure is given in | Appendix A| there are some subtleties 
involved in the choice of grid and interpolation method. An example of the numerical solution for / = 0.1 is 
presented in Figures 1(a) and 1(b) The power- law behaviour at large \z\, with the exponent given by ©, 



is confirmed numerically for / = 0.1 (right panel of Figure 1(a) ) and a range of other values. The iterations 
converge in the negative half-plane, Re(z) < 0, despite the complicated structure of g(z), illustrated in 
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(a) linear scale 



(b) log-linear scale 



Figure 2: The steady-state wealth probability distribution function p s (w) obtained by inverting the Laplace transform g(z) for 
the following values of the transfer fraction: / = 0.5 (bold solid curve), 0.25 (dash-dot curve), 0.1 (dotted curve), 0.05 (dashed 
curve), and 0.025 (thin solid curve). 



Figure 1(b) as it gradually approaches e 
the initial function go{z); e.g. go(z) = e~" 



z for decreasing values of /. The convergence does not depend on 
works just as well for small values of /. 



4. Steady-state wealth distribution by Laplace inversion 

The steady-state probability distribution function p s (w) can be obtained by inverting its Laplace trans- 
form g(z) numerically. A number of inversion algorithms were reviewed recently in [14[ and jl5j]. The 
reviewers advised that at least two different algorithms should be used as a cross-check, because different 
algorithms work well for specific classes of functions and none of the algorithms is universally accurate. 
Fortunately, the algorithms are easy to implement. We test four (referred to as the Euler, Talbot, Stehfcst, 
and Zakian algorithms in 14], [15() and find that the first two give accurate results over a wider range of w. 
Euler has an additional advantage over Talbot: it samples g(z) in the positive half-plane only, where the 
function g(z) has a simpler structure, as one sees in Figu re |l(b)| The results of the inversion are presented 
in Figures [2(a)] and |2(b)| for 0.025 < f < 0.5 and Figures [3(aj| and |3 (b) | for 0.5 </ < 0.9. The exponential 
analytic solution is recovered numerically for / = 0.5. From the output of the inversion algorithms, we 
compute the moments of the distribution ({j,q, /ii, and \xi) and find agreement with the analytical results to 
8 significant digits. In Figures 4(a) and |4(b)] for / = 0.95 and / = 0.05 respectively, we compare the wealth 
distributions obtained from the Laplace transform (curves) and from the agent-based simulations (crosses). 
We find excellent agreement between these two methods for all values of / that we consider. 

The algorithms that perform the inverse Laplace transform suffer from truncation errors. Maximum 
precision is achieved near the peak of p s (w). In the tail, the precision decreases until the results are 
completely dominated by the truncation errors below a threshold value of p s (w). For example, we perform 
all computations with 16 significant digits and achieve ~ 8 significant digits of precision at the peak of 
p s (w), but the algorithms break down at p s (w) < 10~ 8 . 
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(a) log-log scale (b) log-linear scale 

Figure 3: The steady-state wealth probability distribution function p s (w) obtained by inverting the Laplace transform g(z) for 
the following values of the transfer fraction: / = 0.5 (bold solid curve), 0.6 (dash-dot curve), 0.7 (dotted curve), 0.8 (dashed 
curve), and 0.9 (thin solid curve). 

The wealth distributions that we find for / < 1/2 are characterized by power-law behaviour at w <C 1, 
in accord with the analytical results, and approximately exponential tails at large w. A careful examination 
of the tails confirms that the asymptotic behaviour at large wealths is not exactly exponential. However, we 
have not been able to find a closed- form expression for it. The distribution becomes tightly concentrated 
around its peak as / decreases; the peak of the distribution gradually shifts towards w = 1. On the other 
hand, the peak shifts towards w = as / increases; the distributions eventually turns into an exponential 
function for / = 1/2. This overall behaviour is similar to that observed in the asset exchange models based 
on the saving propensity with < A < 1 Q. 

The structure of the steady-state solutions for / > 1/2 is very different. The asymptotic approximation 
p s {w) oc ui a-1 , with a = — l/log 2 (l — /), is valid for / > 1/2 as well and indicates that the distribution 
diverges at w = (as a < 1). For values of / sufficiently close to 1, the wealth distribution acquires a 
shape that is akin to a power law, p s (w) oc w , with overlaid oscillations that become more prominent as 
/ increases. This power-law behaviour cuts off exponentially at some critical wealth that increases slowly 
as / approaches 1. At the same time the exponential drop-off at large wealths becomes shallower as evident 
from Figure [3(b)| The oscillations of p s (w) appear to be periodic on a logarithmic scale, with the period 
depending on /. For example, the periods for / = 0.9 and / = 0.99 are roughly one and two decades 
respectively. This is directly related to the fact that all givers retain 1% of their wealth for / = 0.99 and 
10% for / = 0.9. 

5. Discussion 

We now use the Laplace transform tools developed in Section @] to address two questions that are costly 
to explore with agent-based simulations: the nature of disorder (entropy) and its evolution in the giver 
scheme, and the degree of inequality in the steady state. 
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(a) / = 0.95 (b) / = 0.05 

Figure 4: The population distribution n(w), shown with crosses, as a function of wealth to, measured in ficticious monetary 
units (m.u.) used in the agent-based simulations. The distribution is computed as the number of agents in every unit wealth 
interval after 100 steps in the simulation of the giver scheme with total number of agents N = 4 X 10 5 and transfer parameter 
[(a)] / = 0.95 and [(b)]/ = 0.05. The initial distribution is uniform in the wealth interval [(a)] [0, 100 m.u.] and |(b)| [0, 500 m.u.] . 
The corresponding solution of the steady-state master equation for the same / is shown with a solid curve, with p s (w) scaled 
to conform with the definition of n(w) according to Np s (w/ (w))/ (w) where (w) is the mean wealth. Both the agent-based 
simulations and the master equation predict oscillations in the wealth distribution in | (a) | but not in |(b)| 

5.1. Entropy and the approach to equilibrium 

According to Boltzmann, states with higher entropy are more probable because they correspond to a 
larger number of microscopic configurations of the system. A closed system evolves to a state of maximum 
entropy, i.e. maximum disorder, which is characterized by the Boltzmann-Gibbs distribution. The exponen- 
tial distribution observed in models where the transfer amount is fixed and constant (see section II. C in [l[) 
is thus consistent with entropy maximization ideas. On the other hand, it is argued in [l[ that multiplicative 
asset exchanges may lead to non-exponential distributions because of the broken time-reversal symmetry, 
whereas, in models with fixed additive exchanges, the time-reversal symmetry is preserved. Despite this, the 
entropy maximization technique has been applied in [l6| to an asset exchange model described by a Hamil- 
tonian quadratic in wealth variables. It predicts a gamma distribution of wealth, but p s (w) in multiplicative 
asset exchange models is not a gamma distribution in general. 

In this section, we explore the applicability of the Boltzmann entropy to the giver scheme. Using the 
Laplace transformed solutions of the master equation (fT]) , we compute the steady-state Boltzmann entropjQ 

/>oo 

S s = - dw p s (w)\og[p s (w)] (11) 
Jq 



4 We stress that this definition applies to a normalized distribution with unit mean. Unlike the case of discrete probability 
distributions, continuous entropy can be negative and it is not invariant with respect to the change of variable. It may be 
more appropriate to consider the Kullback-Leibner divergence D s = f^°dw p B {w) log[p s (w) /m(w)], which is a measure of 
the divergence between p s (u>) and the reference distribution m(w). It is convenient to take m(w) = e~ w , in which case the 
divergence essentially reverts to Boltzmann entropy because D s = 1 — S s - 



7 




Figure 5: |(a)| Boltzmann entropy S s of the steady-state distribution as a function of the variance <j\ = f/(l — /). The critical 
values <t^ = 0.062, 1, and 5.098, corresponding to / = 0.058, 0.5, and 0.836 respectively, are indicated with the dotted lines. |(b)| 
Entropy as a function of time for the initial distribution given by il 1 21) with / = 0.058, computed from the multi-agent simulation 
of the giver scheme. For the simulation, the distribution J 1 2D was scaled up to give N = 337123 agents in < w < 1421. To 
compute the entropy, the population distribution produced by the simulation was normalized to a probability distribution with 
unit mean. 



for several values of / in the range 0.01 < / < 0.9 and plot the results in Figure 5(a) Entropy maximization 
arguments (l~7l | imply that the entropy defined by leads uniquely to the exponential distribution, p(w) = 
e~ w with S s = 1, if the only condition is that the mean of the distribution is fixed to \i\ = 1. In our model, 
however, the transfer fraction / places additional constraints on how the distribution of wealth evolves with 
time. Therefore, it is not surprising that we observe a range of steady-state entropies S a ^ 1 corresponding 
to different values of /. The exponential distribution for / = 1/2 appears to be the most disordered state of 
the system with the highest entropy S a = 1. For all other values of /, the entropy is smaller and it becomes 
negative as / approaches or 1. Reading off the graph, we find that S s is negative for < / < 0.058 and 
0.836 < / < 1. The entropy tends to negative infinity as / — > or / — s- 1, which is in accord with the 
behaviour of p s (w) in these limits. 

Negative values of S s are already a warning that the Boltzmann entropy may not be a faithful measure of 
disorder in a multiplicative asset transfer system like the giver scheme. However, the situation worsens when 
we look at how S(t) = — J dw p(w, t) \og[p(w, t)] evolves with time by conducting multi-agent simulations. 
In many cases, it decreases instead of increasing. For example, if we choose p(w, 0) = e~ w initially, S(t) 
decreases with time for all / ^ 1/2 and remains constant for / = 1/2. Moreover, we can easily find realistic 
situations where S(t) does not change monotonically with time, as the experiment described below shows. 
Consider an initial distribution of the form 

( Pi, 0<w<l, 
p{w,0) = < p 2 , K w < w 2 , (12) 
[ 0, otherwise, 

with the parameters pi, p 2 , and w 2 chosen to give 5(0) = (pi ps 0.296, p 2 ~ 1.669, and w 2 re 1.421). The 
evolution of entropy for this initial distribution with / = 0.058 is plotted in Figure [5 (b)| The entropy grows 
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Figure 6: |(a)| The limiting probability distribution £(u>), shown with crosses, obtained from one realization of the asymmetric 
random process fEfl for / = 0.05 after 10 6 iterations. The mean is 0.9998 and the variance is 0.0519 (cf. 0.0526 theoretically 
from the master equation). The steady-state distribution Ps("w) for the same / obtained by Laplace inversion is shown as a 
solid curve for comparison. |(b)|The first 1000 values of {wi}. 



initially but after about ten steps in the simulation it turns over and begins to decrease, eventually reaching 
S s = as expected for / = 0.058. This is in marked contrast to the behaviour of S(t) in an ideal gas, where 
one has dS(t)/dt > according to Boltzmann's H-theorem. 

The population distribution is determined by dividing the wealth axis into small bins and computing the 
number of agents that fall in each bin. One can define the multiplicity W as the number of permutations 
of the agents between different wealth bins such that the occupation numbers of the bins do not change. 
The definition of entropy, S s = logp^, leads to the expression (fTT]) in the continuous limit. Entropy maxi- 
mization under the assumption that total wealth is conserved gives an exponential distribution. However, 
this ignores the global constraints on the probability distribution p s (w) imposed by the transfer fraction 
/. The maximization procedure must take these constraints into account to derive the steady-state wealth 
distribution appropriate to the giver scheme. Unfortunately, at the time of writing, we have been unable to 
derive these additional constraints, and they do not appear in the literature. 

5.2. Random process 

The evolution of wealth in the giver scheme can be analyzed in terms of a random process defined by 

Wi+i = Wi + AiVi, (13) 

where w± = 1 and Awi = +/ or Awi = —fu>i with equal probability. This process is asymmetric, i.e. 
multiplicative in the negative direction and additive in the positive. To illuminate the relationship between 
the giver scheme and the random process we note that (1) an agent's loss of wealth is always proportional 
to his wealth, i.e. it is multiplicative, and (2) an agent's gain of wealth can originate from any other agent 
in the population and therefore equals f(w) on average, or simply / if we set (w) = 1. We compute the 
limiting distribution £(w) of this process by applying (|13|) a sufficiently large number of times and then 
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constructing a histogram of all {wi}. By the ergodic assumption, this is equivalent to computing a large 
number of realizations of this random process and using the final values in each realization to find the 
limiting distribution. 

Figures 6(a) and |6 (b)] display one particular realization of the random process (|13p and the corresponding 



limiting distribution. Given the close link between the transfer model and the random process, it is not 
surprising that the limiting distribution £(w) of the random process (| 13[) appears to be identical to the 
steady-state distribution p s (w) of the giver scheme. We obtain similar results for other values of / <C 1. 
Note that the slight discrepancy between £(w) and p s (w) at w > 2 is due to insufficient sampling of large 
wealths by the random process. The agreement improves as the number of iterations increases. 

5.3. Inequality of wealth 

A traditional measure of inequality in economic systems is the Gini coefficient, defined as G = 1 — 
2 Jq L(X)alX , where L{X) is the Lorenz curve. For a continuous distribution, we have L(w) = /J" dw' w'p s (w'), 
X(w) = Jq dw' p s (w')H and hence 



/>oo pW 

G s = 1 — 2 / dw p s (w) / dw' w'p s (w'), (14) 
Jo Jo 

such that G s = corresponds to perfect equality, and G 8 = 1 to perfect inequality 

In Figure[7]we plot G s versus the steady-state variance u 1 = //(l — /) for 0.01 < / < 0.9. As expected, G s 
increases monotonically with a 2 , since both quantities are measures of dispersion. Interestingly, however, 
there is an inflection point in the G s (a 2 ) curve at a 2 = 1, G s = 1/2, corresponding to the exponential 
distribution (i.e. / = 1/2). For / — > 0, we have p s (w) — > S(w — 1), which corresponds to perfect equality 
since all agents have the same wealth. On the other hand, for / — > 1, the distribution p s (w) becomes sharply 
peaked near w = 0, while the standard deviation approaches infinity. This corresponds to the situation where 
most agents have zero wealth, except for one who has everything, i.e. perfect inequality. 

We can understand the evolution towards inequality in terms of the state vector of the system. Consider 
N agents whose wealths are characterized by random variables Wi, i = 1,2,. . . ,N. The state of the system 
can be described by the phase-space vector w = {w\,W2 1 ■ ■ ■ , tujv)- The constraints that define the phase- 
space are (1) < Wi < 1 for all i, and (2) J^iLi Wi = 1 (we assume for convenience that the total 
wealth is unity). These constraints define a segment of the (N — l)-dimensional hyperplane embedded 
in A^-dimcnsional space. Without any additional constraints, entropy maximization 17| gives gi(wi) = 
(N — 1)(1 — Wi) N ~ 2 for the probability distribution of the wealth Wi of the i-th agent, with mean (wi) = -h 
and variance o~ 2 0i — jyjfeqrn ■ However, in our system, the asset transfer process and the value of the 
parameter / place additional restrictions on the evolution of w. For the increment of the phase-space vector 
w from time tk to time tfc+i, i.e. after one generation of asset transfers, we have 

N 

||Aw|| 2 = J2b»i(tk+i) - w t {t k )] 2 . (15) 

i=l 

The terms in the sum on the right hand side can be split into two groups, associated with the givers and 
the receivers. Since the transfer amount is proportional to the giver's wealth, we get 

||Aw|| 2 = 2/ 2 J2 M*k)] 2 - (16) 

i£ givers 

Therefore the following inequality must always be satisfied: 

||Aw||<2 1 / 2 /l|w||. (17) 



5 Note that we have L(X) < X for all X because L(0) = 0, L(l) = 1, and dL/dX is a monotonically increasing function. 
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Figure 7: Gini coefficient of the steady-state distribution p a (w) as a function of the variance (rj = //(l — /)• 



In addition, we have 



N-V* < llwll < 1 



(18) 



due to the restrictions of the phase-space itself. Note that the state Wi = 1/N for all i is the nearest point 
to the origin. 

When / is small, the norm of the increment ||Aw|| is also small compared with the maximum linear 
extent of the phase space (which equals 2 1 / 2 ); the evolution of w is gradual. Furthermore, ||Aw|| is also 
constrained by ||w||, which can be very small if N is large and all agents are clustered as close to the origin 
as possible. So, if the dispersion in wealth is modest, w moves slower through the phase space than if there 
is great inequality. On the other hand, w changes more rapidly on average if ||w|| is close to unity, which 
corresponds to large inequality. The states of equality are therefore more probable, which explains why the 
steady-state distribution tends towards a delta function for / — > 0. Even if the initial wealth distribution is 
very unequal, w drifts quickly towards the states of near equality. 

When / is close to 1, || Aw|| can be comparable to the size of the phase space. Since / is large, the gains 
in wealth of the individual agents can be large as well. This leads to the situation where a few agents own 
most of the wealth. These agents retain their large wealth for a short time only (typically a few time steps) 
before they become givers and pass their large wealth to other agents. In the extreme case / = 1, one agent 
possesses all the wealth at any instant, while all the other agents have zero wealth. This maximum wealth 
is passed from agent to agent frequently. This corresponds to w jumping from one corner of the phase space 
to another. For / w 1, w evolves similarly, with p s (w) peaking strongly at zero wealth. 



6. Conclusions 

We develop a new technique for computing the steady-state probability distribution of a multiplicative 
asset transfer model, which we call the giver scheme, by Laplace transforming the associated master equation 
to give a functional equation for the characteristic function of the distribution. In the giver scheme, the 
transfer amount fw g is proportional to the giving agent's wealth w g , so the model depends on a single 
parameter / G (0, 1). We develop an efficient iterative method to solve the functional equation for any /, 
and we employ several Laplace inversion algorithms to recover the steady-state distribution p s (w). 
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We comprehensively explore the dependence of the wealth distribution on the value of /, especially 
the thinly studied regime 1/2 < / < 1. We find a stark qualitative difference between the distributions 
for / w (sharply peaked distribution centred around the mean wealth) and / » 1 (broad distribution 
of approximately power-law shape with overlaid oscillations). These two extremes correspond to near- 
perfect equality and inequality respectively, as characterized by the Gini coefficient. Both extremes are also 
characterized by negative Boltzmann entropy. While the regime / is generally thought to represent to 
some extent the exchange processes occuring in the real economy, the regime / ps 1 is probably less applicable 
to realistic economic systems, except perhaps in situations involving extreme leverage. The regime / ~ 1 
may also be relevant to the analysis of gambling, where transitory fortunes are made and lost frequently 

We show that the Boltzmann entropy is unlikely to be a faithful measure of disorder in a multiplicative 
asset transfer system, since it does not vary monotonically as a function of time, assuming the second law of 
thermodynamics. This is an important and counterintuitive result, because the system in the giver scheme 
is closed and the microscopic transfer rules conserve wealth, in a manner reminiscent of the microcanonical 
ensemble in statistical mechanics. In a multiplicative transfer system, the correlations between various 
subsystems (e.g. subclasses corresponding to a particular historical sequence of giving and receiving) and 
the time-reversal asymmetry of the microscopic rules are crucial to the system's dynamics and, therefore, 
cannot be ignored. 



Appendix A. Iterative procedure 

We assume that the computations are carried out with 16 significant digits. For a given complex argument 
z, define a uniform grid u = {uk\i that covers the interval [— 4, log 10 (|z|)]. The approximation ^ gives 
sufficient precision for \z\ < 10 -4 for computations with 16 significant digits. Choose the number of points 
K such that there are a large number of points in every unit interval, say, 10 3 logarithmic grid points per 
decade in [10 , \z\\. Define two auxiliary grids, uro = log 10 (/) + u and = log 10 (l — /) + u, and 

initialize the iterations with gQ{10 u z/\z\) = l/(l + 10"z/|z|). For a given set of values gi(lO u z/\z\), defined on 
the grid u, find the corresponding values on the auxiliary grids by performing a spline interpolation or using 
the approximation ([5]) where appropriate. Then use these values in equation (|10p to find gi+i(10 u z/\z\). 
Continue iterating until the convergence criterion is met (we find that the convergence spreads gradually 
from zero to \z\). Typically the convergence requires a few dozens of iterations for \z\ ~ 100 in the positive 
half-plane. Once the convergence is reached, apply a spline interpolation to find g{z') for any z' along the 
same direction in the complex plane as z, provided that \z\ < \z\. 

The obvious disadvantage of the procedure outlined above is that it relies on interpolation. Its precision 
is therefore limited by the number of points in the grid u, i.e. the discretization of the interval [0, \z\]. It is 
possible, however, to avoid interpolation altogether by defining a special non-uniform grid that is invariant 
with respect to division by / and (1 — /). This gives rise to an alternative procedure for computing the 
iterations. 

Define a grid r k , m = f k (l - f) m with < k < K and < m < M, where K = [log(10" 4 /|z|)/ log(/)] 
and M = [log(10- 4 /|z|)/log(l - /)] are defined such that \z\r Kfi < 10~ 4 and \z\r , M < 10~ 4 . The 
function g{z) on the grid Zk }T n = rk tm z, defined according to gk, m = g(rk,mZ), has the following properties: 
g{fzk,m) = g(Zk+i,m) = 9k+\,m and g[(l - f)z k>m ] = g(z k>m+ i) = gk,m+i- Therefore the iteration rule 
becomes 

9k,m+l / it -.x 

9k,m = ~ , (A.l) 

1 ~ 9k+l,m 

for < k < K — 1 and < m < M — 1. For gx,M one can use the approximation ©. In fact, the Taylor 
expansion can be used for any point |zfc. m | < 10~ 4 . Thus, no interpolation is required and the iterations can 
be computed more efficiently. However, unlike the approach based on interpolation, this procedure must be 
repeated for different arguments even if they lie in the same direction in the complex plane. 
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